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INTERACTION BETWEEN A CIRCULAR INCLUSION 
AND AN ARBITRARILY ORIENTED CRACK* 

by 

F. Erdogan , G. D. Gupta, M . Ratwani 
Lehigh University, Bethlehem, Pa. 18015 


Abstract 

The plane interaction problem for a circular elastic inclusion 
imbedded into an elastic matrix which contains an arbitrarily 
oriented crack is considered. Using the existing solutions for 
the edge dislocations [6] as Green's functions, first the general 
problem of a through crack in the form of an arbitrary smooth arc 
located in the matrix in the vicinity of the inclusion is formu- 
lated. The integral equations for the line crack are then obtained 
as a system of singular integral equations with simple Cauchy 
kernels. The singular behavior of the stresses around the crack 
tips is examined and the expressions for the stress intensity fac- 
tors representing the strength of the stress singularities are 
obtained in terms of the asymptotic values of the density functions 
of the integral equations. The problem is solved for various 
typical crack orientations and the corresponding stress intensity 
factors are given. 


1 . INTRODUCTION 

In fracture studies of ceramics and other composite materials 
it is generally conjectured that the fracture of the solid will 
initiate at and will propagate from a "dominant flaw". This may 
be a manufacturing flaw, it may be caused by residual stresses or 
some other type of loading before the part is put into use, or it 
may result from the growth of a "micro flaw" due to cyclic nature 
of the operating stresses. In some cases it may be possible to 
detect such flaws by using nondestructive testing techniques. 


★ 
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V More often, in studies relating to structural i ntegri ty and rel i - 

ability one simply assumes their existence. Thus, preliminary to 

the application of the relevant fracture theory, one needs to 

solve the mechanics problem for the composite medium consisting of 

inclusions and/or pores and the surrounding elastic matrix which 

i contains the flaw. Generally the flaw is assumed to be an internal 

| ' . • . ... . * ■ 

crack the size of which is of the same order of magnitude as that 

of the inclusions. The exact elasticity treatment of the three 

dimensional problem with a regular or a random array of elastic 

inclusions imbedded into an elastic matrix containing an internal 

crack appears to be hopelessly complicated. However, one may have 

some idea about the response of the composite solid if the solution 

of the related two-dimensional elasticity problem were to be 

avai 1 abl e . 

In this paper such an idealized problem will be considered. 

It will be assumed that the elastic matrix contains only "sparsely" 
distributed inclusions. Hence the mechanical interaction is pri- 
marily between an isolated inclusion or a hole and a line crack 
arbitrarily located in the neighborhood of the inclusion in the 
surrounding elastic matrix which is assumed to be infinitely large. 
The body will be assumed to be in a state of plane strain or gener- 
alized plane stress. The special cases of this problem in which 
the geometry of the medium and the external loads contain a plane 
of symmetry were considered in [1-4]. [4] also includes the solu- 

tion of the problems for two collinear cracks (one in the matrix, 
one in the inclusion), for a crack terminating at the inclusion- 
matrix interface, and for a crack crossing the interface. 
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2. THE INTEGRAL EQUATIONS 


As in many crack problems, the solution of the present problem 
too may be obtained through the superposition of two solutions. 

The first refers to the simple problem of a circular elastic inclu- 
sion inserted into a matrix without the crack. This problem is 
solved under the given system of external loads. In the second 
problem only the stress disturbance due to the existence of the 
crack in the matrix is considered. In this problem the only exter- 
nal loads are the crack surface tractions which are equal in magni- 
tude and opposite in sign to the stresses obtained in the first 
problem along the line which is the presumed location of the crack. 
The nonhomogeneous medium may be subjected to an arbitrary set of 
external loads (including quasi-static thermal loads) applied to 
the matrix and the inclusion. However, it is assumed that the 
dimensions of the matrix are sufficiently large so that in the 
second problem the interaction between the outer boundary of the 
matrix and the crack-inclusion combination may be neglected. "Thus, 
in the second problem, which contains the singular part of the 
solution, the matrix will be assumed as being i nf i ni te . It wi 11 
henceforth be assumed that the solution of the first 'problem is 
known. For example, if the matrix is infinitely large and is sub- 
jected to biaxial stresses at infinity, this solution may be 
obtained by adding the results for two perpendi cul ar uniaxial 


stresses given in [5], namely 
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rr 2 u 2 u 
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(1 + 3^ + 36Rl)s in2 e , 
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2 ( y 2 - y ] ) 
y l + K 1 y 2 * 
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2y 2 + y 1 ( Kg- 1 ) 


y 2~ y l 
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(2. a-c) 


where a Q is the uniaxial stress at infinity, R is the radius of 
the inclusion, r and 8 are the polar coordinates (r is measured 
from the center of the inclusion, 0 is measured from the direction 
of loading, a Q ) , and p-j , k-j and are, respectively, the elastic 

constants of the matrix and the inclusion. 


The integral equations for the stress disturbance problem may 
be obtained by using the solution for the edge dislocations given 
in [6] as the Green's functions. Referring to Figure 1(a), let 
the matrix 1 contain two edge dislocations at the point (x=c o ,y=0) 
with Burgers vectors b v and b u . Following [6], the stresses at a 
point P(x,y) in the matrix may be expressed as 

IT ( K-. +1 ) I 

r °ww(x.y>cj = h (x ,y ,c ) b + h -(x,y,c )b, , 


’o /u x yy 2 


0 ' y 


^(k-j+I ) -| 


°xx (x > y ’ c o> =■ h X xl( x - y - c o )b x + b xx2 (x - y - c o )b y > 


■*! ( K 1 + 1 ) I 

~ °xy (x * y ’ c o ) = h xyl (x - y - c o )b x + h xy2 (x ’ y ’ c o )b y ’ 


(3.a-c) 


where 
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1 +IT1K- 


m( 1 + k i ) 


k 2 - m<i 
K 2 + m 


m = 


X 1 = x “ c o 


x 2 = X 


2 2 2 
r - x £ + y 


r 1 2 = ( x " c 0 ) 2 + y 2 ’ 


R 2 2 ? 

= (x -§-) + y 2 

c o 


Let us now assume that the matrix contains a cut along a smooth 
arc L going through the point (c Q ,0) (Figure la). Let o n (s) and 
a t (s), (s6L) be the normal and tangential components of the stress 
vector on L in the composite without the cut and subjected to the 
given set of external loads, where s is the arc length measured on 

i 

L (i.e., the solution of the first problem). Let the point (c Q ,0) 

correspond to s = s q on L. Once L is specified, from (3) we may 

also obtain the normal and tangential stresses on L due to the^— 

dislocations b and b only as follows: 

^ y 

a n ( s , s_ ) = h , (s >s )b + h 9 (s ,s. ) b ■ , 


o' x 


o' y 


a t (s ’ s o ) = h tl (s ’ s o )b x + h t2 (s ’ s o )b y • 


(7 . a ,b) 


where h . and h t . , ( i =1 , 2 ) may be obtained from (4) by using the 
coordinates on L. If n and t refer to the positive normal and 
tangent to L at s (Figure la) and if a = a(s) is the angle between 
the x-axis and the normal n, noting that (3) is also valid for 
(x,y) € L, o b and a b are given by 

o b = cos 2 a + sin 2 a + 2a| sinacosa , 

n xx y y xy 


b 1 1 12 2 

c° = (Oy y - a^ x )sinacosa + a xy ( cos a - sin a) . 


(8 . a ,b) 
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Substituting x=x(s), y = y(s) from the equation of L, (3), (4), 

(7) and (8) give h ni - and h ^. , ( i = 1 , 2 ) . 

If v/e now assume that b and b v are continuously distributed 

x y 

(unknown) functions (of s Q ) on L, the stress disturbance problem 
(the second problem) may be formulated as 

-o n (s) = ; Eh nl (s,s 0 )b x (s 0 ) + h n2 (s.s 0 )b y (s 0 )]ds 0 -, 

-cr t ( s ) = / [h tl ( s » s o ) b x ( s o ) + h t2^ s,s o^ b y^ s o^ ds o * (9*a,b) 

The si ngl e-val uedness condition of displacement vector requires 
that the density functions of the system of integral equations (9) 
satisfy the following relations: 

I b x ( s 0 ) ds o = 0 > / b y (s 0 ) d s 0 = 0 . ( 1 0. a , b ) 

Superimposed on the solution of the corresponding "first problem", 
in principle (9) and (10) give the solution to any problem of an 
arbitrarily oriented arc-shaped crack in the matrix. 


In particular if the crack is along a straight arc L shown 
in Figure lb the transformations would be considerably simplified. 
Let (t,w) be the fixed coordinate system and the end points of the 
cut be at (t-|,c) and ( 1 2 > c ) (Figure lb). Note that on L 


. b b 

a = a . t = -s , a,,,, = cr = a , cr. = -cr = -a, . 

0 ww nn n wt ns t 


( 11 ) 


The relations and quantities necessary to express the kernels in 
(9) explicitly in terms of t and t Q may then be written as 


sma = 


cosa = 




C c + t, 


^77 


x = t sina + c cosa, y = c sina - t cosa , 
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ww 








yt -t 0 ) 
^ 2 


tt + c 2 - R 2 
0 

42T7T 


r 2 = c 2 + t 2 , 


^ = (t 0 -t ) 2 , 


2 (tt + c 2 - R 2 ) 2 + c 2 (t 0 - t) 2 
r 2 Z = 5 2 £ . (1 2) 

c +t o 


If we now let the stresses on L obtained from the solution of 
the problem without the crack be 


“ ff ww (t ’ c) = P^t) , 


a wt( t,c ) " 


(13) 


and define 


b x' s o> = - f l (t o> • 


b y (s o> = - f 2 (t o> • 


(14) 


the system of integral equations (9) and the conditions (10) may 
be expressed as 

+1) l l 2 

-2.7, Pi' 1 ’ * { J, K ij' t * t o) f j' t o) dt o • . 

z 2 3 ” ’ ' 

(i=l ,2; t 2 <t<t 1 ). (15) 


f 


/ f • (t)dt = 0 , (1=1,2) . 

t 2 


(16) 


It is not difficult to show that at t=t Q the kernels Kjj(t,t 0 ) 

have Cauchy type singularities. In fact, examining the terms 

2 2 

containing x-j/r-j and y/r-j in (4) we find 


.) ! .J— = (. 

' t. - t ’ „ 2 1 


... " . .J t-’rt • ' ,7 > 

'c 2 + t 0 2 


whereas all the remaining terms can be shown to be bounded in the 
closed interval t 2 <_ ( t ,t Q ) £ t-j . Separating the singular parts of 
the kernels, (15) may be written as 


7r( K, + 1 ) 

—2jr — p^t) = / ■ f 1 (t 0 ) 


4: 2 + t, 


dt U 1 

== t 0 - t + l k n (tst o^ f l <V dt o 


+ 1 W t=; tT^t + i ■'k 1 2(t.t 0 )f 2 ( t o)atb •' 


i 2 ' o ' ,- 9 - — s- t - t ‘ i 

t 2 /c 2 + t n 2 0 


tt(ic, + V) . r l 

2y p 2 (t ) = ( MV 


t n dt C 1 

- zp — ■ t~~t + / V^’VMV'M : 

'c 2 + t 2 0 *2 


“ { f 2 ^ 1 0 ^ J-J 2 t n 

*2 If /c 2 + t 0 2 0 


dt n 1 

j— j- ♦./: k 22 (t,t 0 )f 2 (t 0 )dt 0 , 


(t2<t<tj ) , (18. a ,b) 


where k ^ ^ ( t , t Q ) s (i,j=l,2) are bounded functions in the closed 


interval ( t£ ^ { t Q , t) _< t-j ) . For the homogeneous medium u-j=]i 2 » 
kt-| = < 2 » m = 1 » A = 0 = 8 » M = .1/(1 +IC-J ) . Consequently, k.j may be sh 
to vanish. Thus, defining t ,w-components of the dislocations by 


shown 


b x sina - b y cosa = b t = -g ] (t Q ) , 
b x cosa + b y sina = b w = -g 2 (t 0 ) » 
the integral equations (18) may be expressed as 


(19. a, b) 


+1) „ m . g 2 (t O> 

2u, P 1 (t) ~ l t - t dt ( 


t 2 0 


*(<1 + 1 ) 9i (t ) 

2p i P2(t) ' l 2 dt ° 


(t T <t<t 2 ). (20. a, b) 


(20) are the well-known integral equations for a straight crack in 
an infinite homogeneous medium in which the normal and tangential 
crack surface tractions p-| and P 2 are the only external loads and 
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the density functions g-| and g 2 are related to the crack surface 
displacements by [7] 

9 1 ^ = Ft ^ u t " u t^ * ^ 2 ( t (u w - u w ) , (21. a, b) 

where u£ = u t (t,c+0), u” = u t (t,c-0), etc. 

3. STRESS INTENSITY FACTORS 

Following [8] it may be shown that the indexes of the singular 
integral equations (18. a) and (18. b) are +1. Hence, the general 
solution will contain two arbitrary constants which may be deter- 
mined from the sing! e-val uedness conditions (16). Also, the 
fundamental functions of the system are 

w,(t) = w 2 (t) - (t-t 2 )* ,/2 (t,- tr 1/2 , ( 22 ) 

Hence the solution of (18) may be expressed as [8] 

f i (t) = w i (t)F-(t) , ( i = l ,2; t 2 <t<t 1 ), (23) 

where F-j and F 2 are Holder-continuous in the closed interval 
[t 2 ,ti]. In order to investigate the singular behavior of the 
stresses around the end points of the crack, we note that the 
expressions (18) give the normal and the tangential components of 
the stresses in the perturbation problem outside as well as inside 
the cut L , i . e . , 

<J* w (t,c) = Pi^)'* a tw^ t,c ^ = » - (24. a, b) 

for all t in the matrix. Let us now define the following section- 
ally holomorphic functions 
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, , 1 r 1 a ki (t o ,f i (t o ) 

*kj< z) = T V-1 «, 

to 0 


where 


= 7 I 


1 


a k,i<V F i<V 


,1/2 dt o. 


t 2 (t 0 - 2)C(t 0 - t 1 )(t 0 - t 2 )] 

(k,j = 1,2), 


( 25 ) 


l ll 


= -a 


22 




+ t, 


*12 = a 21 


& 


(26) 


+ t, 


-o • v o 

Following [8] the asymptotic behavior of the Cauchy integral 
(25) may be expressed as 

*kj( z * = 1 a k ,(t 2 )F,(t 2 )(t 1 - t 2 ) 1/2 (z - t 2 ) 1/2 


: - a kj( t l) F j{ t l ) ( t l- t 2 )' 1 / 2 (z - t l)' , '‘ + ♦okj' 1 2 ) • 

, . (27) 

where around t-j and t 2 <f> 0 ^ j is generally bounded or at most has 


,- 1/2 


a singularity which is weaker than that of (z-t-)‘ 

0 

<{>kj(z) is holomorphic outside the cut we may write 

1 ^ a ki(^*o^^ i 

U ■ J t - t - dt o - ♦«(*) 

l 2 o 


■ 1/2 


Si nee 


a kj( t 2 )F j (t 2^ 


(t l ■ t 2 ) 1/2 (t 2 - t) 1/2 


a kj^l )F j (t l ^ 


(t r - t 9 ) 1/2 (t- t,) 


+ ^ on ( t ) » 


1/2 v ok j 


(t<t 2 , t>t 1 ). (28) 

From (18) and (28) it is seen that around the crack tips 
t = t -j and t=t 2 the stresses p-j and p 2 will have the conventional 
square-root singularity. If we now define the "stress intensity 
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factors" as follows: 


ki ( t-, ) = 1 im /2( t - t, ) 
t-*-t 1 ' 

P 1 ft) , 



k 9 ( 1 1 ) = 1 i m V2 ( t - t , ) 
t-vt 1 1 

P 2 (t) » 



k-, (to) = 1 im /2(t 9 - t) 
t * >t 2 

P] (t) , 



k ? (t 9 ) = lim /2(t 9 -t) 
t->t 2 c . 

P 2 ( t ) » 


(29. a-d) 

and let 




a = (t 2 - t-, )/2 , 



(30) 

noting that 




/ k 1J (t.t 0 )f j (t 0 )dt 0 ■ 

finite , 

( — 00<t<°°) , 

(31) 


we f i nd 


2m, 


M*! * = " 1 + K 


1 /a(c 2 + t 1 2 ) 


2Pi 


k 2 ^ 1 1 ^ = T+~k 


1 /a(c 2 + t-j 2 ) 


CcF-j(t-|) + ti Fg ( ti ) ] 


[ 1 1 F i ( 1 1 ) - c F 2 ( t -j ) ] 


2y 1 


M 1 ^ = 1 + K 


1 AJ7 


+ -t z ■) 


2m, 


k 2 ( 1 2 ^ “ 1 + k 


^ /a (c 2 + t 9 2 ) 


[cF -| (t 2 ) t 2 F 2 ( t 2 ) ] * 


Ct 2 F-| (t 2 ) - cF 2 (t 2 )] , (32. a-d ) 


where F-j and F 2 are the bounded functions defined by (23). It 
should be obersved that if the transformation (19) is made before 
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m 


solving the integral equations, repeating the foregoing analysis 
or directly from (32) it is seen that in terms of the new unknown 


functions g-j and g 2 the stress intensity factors may be expressed 


as 





2u, 

k 2 ( t 2 ) = lim V2 ( t ~ T 9 1 ( t ) ( 33 . a-d) 

t t 2 1 

where g-j and g 2 are the derivatives of the crack surface displace- 
ments (see (21 ) ) . 


4. NUMERICAL RESULTS 

Once the tractions p-j and p 2 are specified, the system of 
singular integral equations may be solved in a straightforward 
manner (see [9] for an effective numerical technique). After 
solving these equations, all the desired field quantities may be 
evaluated by means of definite integrals with appropriate kernels 
and f-j and f 2 as the density functions. Since the main interest 
in this paper is in the fracture of composite materials the numer- 
ical results will be presented only for the stress intensity 
factors defined by (29) and evaluated from (32). The results are 
given for a uniaxial stress o Q at infinity (see the insert in 
Figures 2-11). The calculated results are shown in Tables 1-7 and 
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Figures 2-11. The stress intensity factor ratios, k.., (i,j=l,2) 

' J 

shown in the tables and the figures are normalized with respect to 
a /a which is the stress intensity factor in a uniaxially stressed 
infinite plane containing a crack of length 2 a perpendicular to 
the direction of loading. Thus - 

Mt.) 

k., * — — — . (i.j = 1,2) , (34) 

c o /a 

where k - ( t - ) is defined by (29). The length parameters a, b, c, 

* J 

and R are shown in the insert of Figures 2-11. In the numerical 
analysis it is assumed that the crack is perpendicular to the 
direction of external loading, a . Only two material combinations 
have been considered. The first refers to a circular hole (i.e., 
^2 = 0) and in the second it is assumed that 

(vj^/p-) ) = 23 , k-j = 1 .6 , <2 - 1 .8 (35) 

which roughly corresponds to a metallic inclusion imbedded into an 
epoxy-type matrix. 

Table 1 and Figure 2 show the results for a radial crack 
perpendicular to o . Note that as the crack tip t 2 approaches the 
hole k^ g° e s to infinity and as it approaches the inclusion- 
matrix interface (where > P-j ) 2 9 0es to zero. This problem 

of a crack terminating at and going through the interface was 
extensively studied in [10, 11, and 4]. It should also be noted 
that for ^2 = 0 b/a = 3, k-j-j will be finite (see [4]). 
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Table 1. The stress intensity factors 
for a symmetrically located radial crack 
{R/a = 2, c - 0). 



^2 

= 0 

y 2 = 

2 3 y i 

b/a 

k n 

■ 

k 12 

k ll 

k l 2 

3.0 


-> CO 


+ 0 

3.2 

1 . 41 7 

2.274 

0.827 

0.467 

3.5 

1.290 

1.722 

0.874 

0.671 

4 

1.188 

1.394 

0.918 

0.821 

5 

1.102 

1 .174 

0.957 

0.924 

6 

1.065 

1.099 

0.973 

0.959 

8 

1.033 

1.045 

0.987 

0.982 

00 

1.0 

1.0 

1.0 

1.0 


Table 2. Stress intensity factors for a symmetrically 
located "tangential crack" perpendicular to the load 
(b = 0, R = 2a). 



Vi 

= 0 

y 2 / y ] 

= 23 

c/a 

k n = k i 2 

k 21 = " k 22 

k ll = k l 2 

k 21 = - k 22 

2.1 

0.0626 

- 0.170 

1.201 

0.105 

2.5 

0.260 

- 0.101 

1.196 

0.171 

3 

0.395 

- 0.0644 

1.196 

0.202 

4 

0.574 

- 0.0175 

1.180 

0.180 

6 

0.768 

0.0274 

1.110 

0.109 

10 

0.906 

0.0361 

1.046 

0.0567 

CO 

1.0 

0 

1.0 

0 
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Table 3. Stress intensity factors for a crack perpendicular 
to the external load (R=2a,b = 3a). 





= 0 



y 2 = 

23p 1 

-j 

c/a 

k n 

k 12 

k 21 

k 22 

k Tl 

k l 2 

k 21 j 

k 22 

0 

mm 

oo 

+ 0 

-*• 0 


- 0 

0 

-»• 0 

0.3 

1 

4.267 

0.072 

-0.391 

0.784 

0.225 

-0.004 

0.072 

0.5 


3.070 

0.113 

-0.412 

0.792 

0.341 

-0.006 

0.101 

1.0 

BBSS 

1.969 

0.177 

-0.316 

0.817 

0.613 

-0.005 

0.057 

1 ,5 

1.357 

1.552 

0.184 

-0.241 

0.839 

0.763 

0.008 

-0.007 

2 

1.299 

1.337 

0.150 

-0.243 

0.860 

0.845 

0.034 

-0.021 

3 

1.189 

1.083 

DSl 

-0.282 

0.905 

0.953 

0.089 

-0.001 

4 

1.091 

0.957 

BBS 

-0.261 

0.951 

1.014 

0.117 

0.002 

8 

0.955 

0.911 

EBB 

-0.114 

1.020 

IBB 


-0.026 

CO 

1.0 

1.0 

0 

0 

1.0 

IBS 

[ 0 " 
I 
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Table 4. Stress intensity factors for a crack perpendicular 
to the external load (R = 2a, c = a) . 




y 2 

= 0 



^2 = 

23 p -| 

. 

b/a 

k n 

k l 2 

k 21 

k 22 

k 

11 

k l 2 

k 21 

k 22 

/3/2 

2.8 


OO 

2.700 

0.206 

->• - CO 

-0.790 

0.781 

.-»■ 0 
0.483 

0.002 

-> 0 
0.167 

3.0 

1.430 

1.970 

0.177 

-0.316 

0.817 

0.613 

-0.005 

0.057 

3.5 

1.274 

1.548 

0.134 

-0.047 

0.878 

0.752 

-0.012 

-0.047 

4 

1.193 

1.364 

0.107 

0.023 

0.914 

0.833 

-0.012 

-0.068 

6, 

1.070 

1.109 

0.048 

0.032 

0.970 

0.952 

-0.006 

-0.041 

8 

1.036 

1.049 

0.025 

0.015 

0.985 

0.980 

-0.002 

-0.021 

10 

1 .022 

1.028 

0.014 

0.007 

0.991 

0.989 

-0.0004 

-0.012 

CO 

1.0 

1.0 

0 

0 

1.0 

1.0 

0 
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Table 5. Stress intensity factors for a crack perpendicular 
to the external load (c = R, b-a = 0.2R; a = constant, R variable, 
or R=constant, a variable). 




^2 

= 0 

' 


y 2 = 

2 3 y i 


R/a 

k n 

k 12 

k 21 

k 22 

k ll 

k l 2 

k 21 

■. 

k 22 

5 

0.366 

0.0582 

-0.102 

-0.0210 

1.120 

1 .096 

0.0256 

0.124 

3 

0.882 

0.257 

-0.0774 

-0.349 

0.984 

1 .063 

0.139 

0.150 

2 

1.247 

0.600 

0.0935 

-0.826 

0.868 

1.022 

0.172 

0.162 

1.5 

1.354 

0.926 

0.246 

-1.226 

0.829 

0.990 

0.168 

0.162 

1.0 

1.312 

1 .384 

0.400 

-1.719 

0.839 

0.951 

0.165 

0.153 

0.75 

1.225 

1.679 

0.453 

-1.972 

0.866 

0.929 

0.176 

0.144 

0.5 

1.092 

2.180 

0.486 

-2.239 

0.899 

0.909 

0.194 

0.129 


Table 6. Stress intensity factors for a crack perpendicular 
totheexternalload(R=2a,c=2.2a). 
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u 2 

= 0 



^2 = 

23^ 


b/ a 

k n 

k l 2 

k 21 

k 22 

k i i 

1/ 

K 1 2 

k 21 

k 22 

0 

0.146 

0.146 


0.113 

1.199 

1.199 

0.127 

-0.127 

0.5 

0.549 

0.008 

1 

-0.101 

1.096 

1.214 

0.202 

-0.026 

1.0 

0.952 

0.196 


-0.372 

0.972 

1.156 

0.212 

0.031 

2 

1.315 

0.926 

0.092 

-0.526 

0.847 

0.994 

0.122 

0.035 

4 

1 .188 

1.268 

0.123 

-0.060 

0.911 

0.871 

0.011 

-0.066 

6 

1.084 

1.126 

0.0765 

0.025 

0.961 

0.941 

-0.0016 

-0.060 
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Table 7. Stress intensity factors for a crack 
perpendicular to the external load (c=R/2). 


P 2 = 0 

R/a 

b/a 

k ll 

k 12 

k 21 

k 22 

■ 

5 

5.5 


2.426 

8 . 22 x 1 0 “^ 

mm 

3 

3.7 

1.721 

2.668 

0.148 

H| 

2 

2.8 

1.562 

2.862 

0.206 

wBm 

1.5 

2.35 

1.448 

2.993 

0.217 

-0.645 

1.0 

1.9 

1.309 

3.139 

0.205 

-0.472 

0.75 

1 .675 

1.232 

3.185 

0.188 

-0.319 

in 

• 

o 

1.45 

1.153 

3.120 

0.165 

-0.0976 

V~2 ~ 23^] 

5 

5.5 


0.508 


HjajgjU 

3 

3.7 

0.716 

0.483 

! : . 


2 

2.8 

0.781 

0.469 


0.156 

1.5 

2.35 

0.827 

0.463 


0.139 

1.0 

1.9 

0.881 

0.460 

0.026 

0.113 

0.75 

1.675 

0.911 

0.464 

0.040 

wmm 

0.5 

1.45 

0.940 

0.481 

0.056 

BISS 














Table 2 and Figure 3 show the results for the other symmetric 
crack geometry, namely, b = 0 and c variable. In this case it is 
seen that the stress intensity factor ratio for the cleavage 
mode may be greater than 1 for the stiffer inclusion and may be 
considerably less than 1 for the hole. This and the similar 
results observed in some of the subsequent figures may at first 
appear to be somewhat paradoxical. However, if one considers the 
distribution of the stresses obtained from (1) in the absence of 
the crack which is shown in Figures 4 and 5, the explanation for 
these trends would be clear. The stresses given in Figures 4 and 
5 are related to the input functions in (18) by 

P-J ( b ) = - a yy(t»c) » p 2 (t) = “ a xy^ t,C ^ ' (36) 

Other results which may be of interest in the applications 
are given in Tables 3-7 and Figures 6-11. Since the shear compo- 
nents, kgj* ( j = 1 , 2 ) of the stress intensity factor are not zero', 
in a brittle or quasi -bri ttl e matrix the crack propagation would 
not be expected to be in the plane of the crack. From the results 
given in this paper, it is not difficult to show that, generally 
for the crack tip near the interface, the crack would propagate 
towards the interface if ^ = 0 or P 2 < ^l ’ anc * away from it if 
P 2 > u-j . A quantitative model for this phenomenon was discussed 
in [12]. 
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Figure 3. The stress intensity factors for a symmetrically 
located "tangential crack" perpendicular to the load (b=0, 
R = 2a ' 
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Figure 8. Stress intensity factors for a crack in the matrix 
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